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(57) ABSTRACT 

A method for registering a first image to a second image using 
a similarity transformation. The each image includes a plu- 
rality of pixels. The first image pixels are mapped to a set of 
first image coordinates and the second image pixels are 
mapped to a set of second image coordinates. The first image 
coordinates of two reference points in the first image are 
determined. The second image coordinates of these reference 
points in the second image are determined. A Cartesian trans- 
lation of the set of second image coordinates is performed 
such that the second image coordinates of the first reference 
point match its first image coordinates. A similarity transfor- 
mation of the translated set of second image coordinates is 
performed. This transformation scales and rotates the second 
image coordinates about the first reference point such that the 
second image coordinates of the second reference point 
match its first image coordinates. 

17 Claims, 6 Drawing Sheets 
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GEOMETRIC REGISTRATION OF IMAGES 
BY SIMILARITY TRANSFORMATION USING 
TWO REFERENCE POINTS 


GOVERNMENT FUNDING 5 


The U.S. Government has a paid-up license in the present 
invention and the right in limited circumstances to require the 
patent owner to license others on reasonable terms as pro- 
vided for by contract as awarded by National Aeronautics and 1 0 
Space Administration under funding numbers NASA Space 
Grant (NNG05G092H) and NASA (NAG5-11773), and by 
the office of Naval Research, under funding number ONR 
(N00014-03-1-0337). 


FIELD OF THE INVENTION 


The present invention concerns a method of registering 
images using reference points. In particular, this method is 
directed to improving registration of images that are partially 20 
occluded, such as satellite images of predominately oceanic 
regions. 

BACKGROUND OF THE INVENTION 

25 

In the case of images with large geometric distortions, an 
orthorectification method may be necessary to reduce the 
distortion due to the relief or the viewing angle of the images. 

For a first order polynomial image registration model, a mini- 
mum of three reference points are required to calculate the 30 
model. For satellite imagery these reference point are known 
as ground control points (GCPs). In practice, it is usually 
better to use more than three reference points. However, it 
may not always be practical to obtain more than three refer- 
ence points, for example, registration of satellite images of 35 
coastal and oceanic regions. 

Present geometric registration methods use an affine trans- 
formation that requires more than three reference points. To 
determine the correct geometric registration of high spatial 
resolution images measured by multi -satellite sensors over 40 
the ocean, it is often desirable to use only images that include 
coastal regions. Still, even these methods may not be appli- 
cable for situations where there are few GCPs available, such 
as in oceanographic studies using satellite multi-sensor mea- 
surements, which often include significant cloud contamina- 45 
tion. The difficulty of registering images is increased by the 
cloud contamination, as the presence of cloud cover that may 
occlude some of the GCPs in one or both of the images. 
Accordingly, it is desirable for improved registration methods 
aimed at studying images of oceanic regions, as well as land 50 
regions. 

The present invention addresses this need among others. 


translated set of second image coordinates is performed. This 
transformation scales and rotates the second image coordi- 
nates about the first reference point such that the second 
image coordinates of the second reference point match its first 
image coordinates. 

Another exemplary embodiment of the present invention is 
a method for registering a first image to a second image using 
a similarity transformation. The each image includes a plu- 
rality of pixels. The first image pixels are mapped to a set of 
first image coordinates and the second image pixels are 
mapped to a set of second image coordinates. A first plurality 
of potential reference points are determined in the first image 
and a second plurality of potential reference points are deter- 
mined in the second image. A set of reference points common 
to the two pluralities of potential reference points is identi- 
fied. The first image coordinates and the second image coor- 
dinates of each reference point in the set of identified common 
reference points is determined. A Cartesian translation of the 
set of second image coordinates is performed such that the 
second image coordinates of the first reference point match its 
first image coordinates. If the set of reference points includes 
exactly two reference points, a similarity transformation of 
the translated set of second image coordinates is performed to 
scale the second image coordinates and rotate the second 
image coordinates about the first reference point such that the 
second image coordinates of a second reference point of the 
set of reference points match its first image coordinates. If the 
set of reference points includes more than two reference 
points, an affine transformation of the translated set of second 
image coordinates is performed to register the second image 
coordinates of the set of reference points to the first image 
coordinates of the set of reference points. 

A further exemplary embodiment of the present invention 
is a computer readable medium adapted to instruct a general 
purpose computer to register a first image to a second image 
using a similarity transformation. The each image includes a 
plurality of pixels. The first image pixels are mapped to a set 
of first image coordinates and the second image pixels are 
mapped to a set of second image coordinates. The first image 
coordinates of two reference points in the first image are 
determined. The second image coordinates of these reference 
points in the second image are determined. A Cartesian trans- 
lation of the set of second image coordinates is performed 
such that the second image coordinates of the first reference 
point match its first image coordinates. A similarity transfor- 
mation of the translated set of second image coordinates is 
performed. This transformation scales and rotates the second 
image coordinates about the first reference point such that the 
second image coordinates of the second reference point 
match its first image coordinates. 


SUMMARY OF THE INVENTION 

An exemplary embodiment of the present invention is a 
method for registering a first image to a second image using a 
similarity transformation. The each image includes a plurality 
of pixels. The first image pixels are mapped to a set of first 
image coordinates and the second image pixels are mapped to 
a set of second image coordinates. The first image coordinates 
of two reference points in the first image are determined. The 
second image coordinates of these reference points in the 
second image are determined. A Cartesian translation of the 
set of second image coordinates is performed such that the 
second image coordinates of the first reference point match its 
first image coordinates. A similarity transformation of the 


BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is best understood from the following 
detailed description when read in connection with the accom- 
panying drawings. Included in the drawing are the following 
figures: 

FIG. 1 is a flow chart illustrating an exemplary method of 
register two images using two reference points according to 
the present invention. 

FIG. 2 is a simplified map illustrating an exemplary coastal 
region with a number of ground control points according to 
the present invention. 

FIGS. 3 A and 3B illustrate exemplary simulated images of 
portions of the coastal region illustrated in FIG. 2. 
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FIG. 4A illustrates the exemplary simulated images of 
FIGS. 3 A and 3B overlaid during the exemplary registration 
method of FIG. 1. 

FIG. 4B illustrates the exemplary simulated images of 
FIGS . 3 A and 3B registered by the exemplary method of FIG. 5 

1 . 

FIG. 5 is a flow chart illustrating an exemplary method of 
register two images according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION to 

Exemplary embodiments of the present invention include 
methods of registering images with only two common refer- 
ence points using similarity transformations. These exem- 
plary registration methods may also include the use of more 15 
accurate registration algorithms when the images have more 
than two common reference points. These exemplary meth- 
ods may enhance present geometric registration processes 
and may also assist in a wide range of applications in image 
analysis. The images used may include satellite or other high 20 
altitude images of a planetary surface. 

Exemplary embodiments of the present invention may be 
used to obtain a geometric registration using only two refer- 
ence points. These exemplary methods include processes of 
translating, scaling and rotating at least one of the images to 25 
match a particular size and position. Present technology for 
the geometric registration, which is also known as georefer- 
encing, georegistration, georectification, and/or geocorrec- 
tion in the geographic and planetary sciences, requires more 
than three reference points and is typically accomplished 30 
using an affine transformation. One advantage of exemplary 
embodiments of the present invention over these methods is 
that exemplary embodiments of the present invention may be 
applicable to wide ranges of imagery, e.g. satellite images of 
offshore regions, which have few potential reference points, 35 
making it difficult to use affine transformation techniques. 
However, exemplary embodiments of the present invention 
may be applicable to these kinds of areas as they require only 
two common references. 

FIG. 1 illustrates an exemplary method for registering one 40 
image to another image using a similarity transformation. The 
images each include a plurality of pixels. The pixels of each 
image may be mapped to a set of image specific coordinates. 
The process of registration involves transforming at least one 
of these sets of image specific coordinates so that the pixels 45 
may be mapped to a common set of coordinates. In the 
descriptions of the exemplary methods of the present inven- 
tion, registration is described in terms of transforming one set 
of image specific coordinates to correspond to the other set 
image specific coordinates. One skilled in the art will under- 50 
stand that this approach has been adopted merely for simplic- 
ity of illustration and is not intended to be limiting. 

The images may illustrate different, overlapping portions 
of a scene. In this situation, registration of the images may 
allow a mosaic of the scene to be formed. Alternatively, the 55 
images may represent data from different sensors, for 
example visible and infrared sensors, that illustrate approxi- 
mately the same scene. Proper registration of such images 
may allow comparison of information about objects in the 
scene detected by the different sensors. 60 

One skilled in the art will understand that exemplary meth- 
ods of the present invention may be used sequentially to 
register many images. It is also contemplated that a combi- 
nation of overlapping images by the same of similar sensors 
and images of approximately the same portion of the scene 65 
taking with different sensor types may all be registered using 
exemplary embodiments of the present invention. 


4 

First image coordinates for two reference points, a first 
reference point and a second reference point, are determined 
in the first image, step 100. Second image coordinates of the 
first and second reference points in the second image are also 
determined, step 102 . It is noted that each reference point may 
correspond to one or more pixels in the images. The origin of 
the set of first image coordinates may be set to one of the 
corners of the first image, the center of the first image, or 
another predetermined location in the first image. Alterna- 
tively, the origin of the first image coordinates may be set to 
the location of the first reference point in the first image. In 
this case, initial first image coordinates of the first reference 
point in the first image are determined. Then a Cartesian 
translation of the set of first image coordinates is performed 
such that the first image coordinates of the first reference 
point are the origin, i.e. (0, 0). The first image coordinates of 
the second reference point in the firstimage are then deter- 
mined using the translated set of first image coordinates. 

One example of imagery that may be used with exemplary 
methods of the present invention high altitude images of a 
planet surface. FIGS. 2, 3 A, 3B, 4A, and 4B illustrate the use 
of the exemplary method of FIG. 1 on simulated high altitude 
satellite images of the Korean peninsula and adjacent seas. 
FIG. 2 shows simplified map 200 of these coastal regions. The 
map includes coastline 202, open water 204 and ground con- 
trol points (GCPs) 206. These GCPs represent easily identi- 
fiable features of the coastline that may be used as reference 
points for registering images. Actual imagery may include 
significantly more GCPs than shown in simplified map 200; 
however, the reduced number of GCPs in FIGS. 2, 3 A, 3B, 
4A, and 4B is used in these figures to help illustrate the 
usefulness of exemplary methods of the present invention and 
to simplify the drawings. This use of a reduced number of 
GCPs, i.e. potential reference points, is not intended to be 
limiting. 

The reference points in the two images may be determined 
by identifying two easily recognizable features that are vis- 
ible in both the first image and the second image. FIG. 3 A 
illustrates first simulated image 300 of a portion of the scene 
illustrated in simplified map 200. FIG. 3B illustrates second 
simulated image 308 of another, overlapping, portion of the 
scene illustrated in simplified map 200. Both of these images 
include some open water regions 204, sections of coastline 
202, and regions obscured by cloud cover 302, as well as 
several GCPs (shown as open circles). Two of these GCPs, 
GCP 304/304' and GCP 306/306', are visible in both first 
simulated image 300 and second simulated image 308. These 
two GCPs may desirably be selected as the first reference 
point, GCP 304/304', and the second reference point, GCP 
306/306'. 

Once the two reference point are selected and the coordi- 
nates in the set of first image coordinates and the set of second 
image coordinates are determined, a Cartesian translation of 
the set of second image coordinates is performed such that the 
second image coordinates of the first reference point match 
the first image coordinates of the first reference point, step 
104 (FIG. 1). FIG. 4A illustrates this step graphically. In 
overlay image 400, second simulated image 308 has been 
translated such that first simulated image 300 and second 
simulated image 308 are overlaid and first GCP 304 from first 
simulated image 300 and first GCP 304' from second simu- 
lated image 308 are coincident, combined first GCP 402. The 
two dashed lines connecting combined first GCP 402 to sec- 
ond GCPs 306 and 306' illustrate differences in rotation and 
scale between simulated image 300 and 308, respectively, 
that are compensated for to complete registration of these 
images. 
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A similarity transformation of the translated set of second 
image coordinates is performed to scale the second image 
coordinates and rotate the second image coordinates about 
the first reference point, step 106 (FIG. 1), such that the 
second image coordinates of the second reference point 
match the first image coordinates of the second reference 
point. 

The first image coordinates of the first reference point and 
the second reference point may be written (x 1 ,y 1 ) and (x 2 ,y 2 ), 
respectively. The second image coordinates of the first refer- 
ence point and the second reference point may be written 
(Xi'yV) and (x 2 ',y 2 ), respectively. A similarity transforma- 
tion may be used to determined scale, A, and rotational angle, 
0, with respect to the first image coordinates to transform the 
second image coordinates to match the coordinates of the 
reference points, i.e. 

(*kPi)“K*i hV) and (x 2 ,y 2 )^(x2,y2) (1) 

Using a rotation matrix representation with a scale, A, and 
a rotational angle, 0, the relation ship of these two reference 
points may be written: 

x 1 '=A cos Qx { +A sin 0jy x and y^-A sin Qx { +A cos 

0^i 

x 2 '=A cos Qx 2 +A sin Qy 2 and y 2 ’=-A sin Qx 2 +A cos 

0^2 (2) 

These equations give the relation between the coordinates 
(x2,y2) of a point relative to the (x2,y2) axes and the coordi- 
nates (x2',y2') of the same point relative to rotated axes (x2’, 
y2'). 

Given that the first image coordinates and the translated 
second image coordinates of both reference points have been 
determined in steps 100 and 104, the scale, A, and the rota- 
tional angle, 0, may be computed from Equation (2). First, to 
simplify these equations, set a=A cos 0 and (3=A sin 0, or: 

A=( <x 2 +|F) 1/2 and O^arT 1 (p/a) (3) 

Using this substitution, Equation (2) may be rewritten in a 
linear form, i.e.: 

andj^ --px^a^ 

x 2 -0LT 2 +pv 2 and y 2 ,= -^x 2 +ay 2 (4) 

It is noted that there are two unknown parameters, a and (3, 
but four equations. Thus, a and (3 are over constrained. How- 
ever a least square method (LSM) may be used to solve the 
system of linear equations. The error, E, may be written by: 

E=((Xx i +$y l -xS) 2 +(-fa i +ay i -ytf+(ax 2 +$y 2 -x^f+ 

(-px 2 +a y 2 -y 2 f (5) 

To determine minimum error to obtain a and (3, i.e. the best 
fit, two partial differential equations are satisfied, i.e., 

3£/Oa=0 and dE^=0 (6) 

and a and (3 are determined to be: 

a=(*i*i '+* 2*2 W (xi 2 +yi 2 +x 2 2 +y 2 2 ) 

P=(-* 1>>1 '+)>!*! -* 2 y 2 , +y 2 x 2 ')/(x 1 2 +y 1 2 +x 2 2 +y 2 2 ) (7) 

These values of a and (3 may be plugged into Equation (3) 
to obtain A and 0. 

Combined image 404 of FIG. 4B illustrates the resulting 
scaling and rotation of simulated second image 308. In this 
image, second GCPs 306 and 306’ are coincident at combined 
second GCP 406. Combined image 404 demonstrates how the 
exemplary image registration method of FIG. 1 may be used 
to form a mosaic from two or more smaller initial images. 
Reduced cloud cover 408 remains only in those portions of 
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combined image 404 that were not visible in either first simu- 
lated image 300 or second simulated image 308. Thus, com- 
bined image 404 is less affected by the cloud contamination 
than either first simulated image 300 or second simulated 
5 image 308. 

It is noted that, in the example use demonstrated in FIGS. 
2, 3A, 3B, 4A, and 4B, two simulated images of relatively 
small portions of the total scene were registered to each other; 
however, it is contemplated that exemplary embodiments of 
to the present invention may also be used to register a relatively 
small image to a reference image that includes a relatively 
large number of predetermined features to use as potential 
reference points. For example, first simulated image 300 
could be registered to simplified map 200. In this case, first 
15 simulated image 300 could be scanned to identify any of 
GCPs 206 to use as reference points. 

FIG. 5 illustrates another exemplary image registration 
method according to the present invention. In this exemplary 
method, the first image includes a first plurality of pixels 
20 mapped to a set of first image coordinates and the second 
image includes a second plurality of pixels mapped to a set of 
second image coordinates. A first plurality of potential refer- 
ence points is determined in the first image and a second 
plurality of potential reference points is determined in the 
25 second image, step 500. Each of the images may include 
imagery of a portion of a scene that includes a set of prese- 
lected features. Determining the pluralities of potential refer- 
ence points may then involve identifying subsets of these 
preselected features of the scene that are visible in the each 
30 image and determining depictions of these subsets of features 
to be the potential reference points. 

Once the pluralities of potential reference points in each 
image are determined, a set of reference points common to 
both pluralities of potential reference points is identified, step 
35 502. This set of reference points may desirably include 
exactly two, three, or four elements. If more than four com- 
mon reference points exist, it may be desirable to only select 
four to avoid having to perform additional calculations; how- 
ever, in some cases it may be desirable to select additional 
40 reference points due to potential distortions that may reduce 
the accuracy of the image registration. 

The first image coordinates and the second image coordi- 
nates of each reference point in the set of reference points 
identified in step 502 are determined, step 504. One of the 
45 reference points is selected to be a first reference point and a 
Cartesian translation of the set of second image coordinates is 
performed such that the second image coordinates of this first 
reference point match its first image coordinates, step 506. 

It is then determined how many common reference points 
50 are in the set of reference points, step 508 . If there are exactly 
two reference points in the set of reference points, a similarity 
transformation of the translated set of second image coordi- 
nates is performed, step 510, to register the second image 
coordinates of the set of reference points to the first image 
55 coordinates of the set of reference points, as described above 
with reference to the exemplary method of FIG. 1 . If there are 
more than two reference points in the set of reference points, 
an affine transformation, of the translated set of second image 
coordinates is performed, step 510, or other prior art trans- 
60 formation method, to register the second image coordinates 
of the set of reference points to the first image coordinates of 
the set of reference points. 

The various exemplary embodiment of the present inven- 
tion may be carried out through the use of a general-purpose 
65 computer system programmed to perform the steps of the 
exemplary methods described above with reference to FIGS. 
1 and 5. Exemplary general-purpose computer systems may 
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include personal computers, work stations, distributed pro- 
cessing computer networks, and parallel processing com- 
puter systems. Parallel or distributed processing may be desir- 
able for substantially real time applications involving the 
substantially concurrent prediction of future quotes for a plu- 5 
rality of financial instruments. Dedicated special-purpose 
computing systems may also be designed for performing 
exemplary methods of the present invention as well. 

Additionally, it is contemplated that the methods previ- 
ously described may be carried out within a general purpose to 
computer system instructed to perform these functions by 
means of a computer-readable medium. Such computer-read- 
able media include; integrated circuits, magnetic and optical 
storage media, as well as audio -frequency, radio frequency, 
and optical carrier waves. 15 

Although the invention is illustrated and described herein 
with reference to specific embodiments, the invention is not 
intended to be limited to the details shown. Rather, various 
modifications may be made in the details within the scope and 
range of equivalents of the claims and without departing from 20 
the invention. 

What is claimed: 

1. A method for registering a first image to a second image 
using a similarity transformation, the first image including a 25 
first plurality of pixels mapped to a set of first image coordi- 
nates and the second image including a second plurality of 
pixels mapped to a set of second image coordinates, the 
method comprising the steps of: 

a) determining the first image coordinates of a plurality of 30 
reference points in the first image; 

b) determining the second image coordinates of the plural- 
ity of reference points in the second image; 

c) performing a Cartesian translation of the set of second 
image coordinates such that the second image coordi- 35 
nates of a first reference point of the plurality of refer- 
ence points match the first image coordinates of the first 
reference point; and 

d) if the plurality of reference points includes exactly two 
reference points, performing a similarity transformation 40 
of the translated set of second image coordinates to scale 
the second image coordinates and rotate the second 
image coordinates about the first reference point such 
that the second image coordinates of a second reference 
point of the plurality of reference points match the first 45 
image coordinates of the second reference point; and 

e) if the plurality of reference points includes more than 
two reference points and no less than two reference 
points, performing an affine transformation of the trans- 
lated set of second image coordinates to register the 50 
second image coordinates of the plurality of reference 
points to the first image coordinates of the plurality of 
reference points. 

2. The method according to claim 1, wherein step (a) 

includes the steps of: 55 

al) determining initial first image coordinates of the first 
reference point in the first image; 

a2) performing a Cartesian translation of the set of first 
image coordinates such that the first image coordinates 
of the first reference point are (0, 0); and 60 

a3) determining the first image coordinates of the second 
reference point in the second image using the translated 
set of first image coordinates. 

3. The method according to claim 1, wherein step (a) 

includes the steps of: 65 

al) identifying a feature visible in both the first image and 
the second image to be the first reference point; 
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a2) identifying another feature visible in both the first 
image and the second image to be the second reference 
point; and 

a3) determining the first image coordinates of the first 
reference point and the second reference point in the first 
image. 

4. The method according to claim 1, wherein: 

the first image is a reference image including a plurality of 
predetermined features; and 
step (a) includes the steps of: 

al) identifying one of the plurality of predetermined fea- 
tures in the second image to be the first reference point; 
a2) identifying another of the plurality of predetermined 
features in the second image to be the second reference 
point; and 

a3) determining the first image coordinates of the first 
reference point and the second reference point in the first 
image. 

5. The method according to claim 1, wherein: 

the first image and the second image are high altitude 
images of a planet surface; and 
the first reference point is a first ground control point 
(GCP) and the second reference point is a second GCR 

6. The method according to claim 5, wherein: 

the high altitude images of the planet surface each include 
a section of coastline; and 

the first GCP is a feature of the section of coastline and the 
second GCP is another feature of the section of coast- 
line. 

7. A method for registering a first image to a second image 
using a similarity transformation, the first image including a 
first plurality of pixels mapped to a set of first image coordi- 
nates and the second image including a second plurality of 
pixels mapped to a set of second image coordinates, the 
method comprising the steps of: 

a) determining a first plurality of reference points in the 
first image; 

b) determining a second plurality of reference points in the 
second image; 

c) identifying a set of reference points common to the first 
plurality of reference points and the second plurality of 
reference points; 

d) determining the first image coordinates and the second 
image coordinates of each reference point in the set of 
reference points identified in step (c); 

e) performing a Cartesian translation of the set of second 
image coordinates such that the second image coordi- 
nates of a first reference point of the set of reference 
points match the first image coordinates of the first ref- 
erence point; 

f) if the set of reference points includes exactly two refer- 
ence points, performing a similarity transformation of 
the translated set of second image coordinates to scale 
the second image coordinates and rotate the second 
image coordinates about the first reference point such 
that the second image coordinates of a second reference 
point of the set of reference points match the first image 
coordinates of the second reference point; and 

g) if the set of reference points includes more than two 
reference points and no less than two reference points, 
performing an affine transformation of the translated set 
of second image coordinates to register the second 
image coordinates of the set of reference points to the 
first image coordinates of the set of reference points. 

8. The method according to claim 7, wherein: 
a set of features in a scene are preselected; 
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the first image includes imagery of a first portion of the 
scene; 

the second image includes imagery of a second portion of 
the scene; 

step (a) includes the steps of: 

al) identifying a first subset of features of the preselected 
set of features of the scene that are visible in the first 
image; and 

a2) determining depictions of the first subset of features in 
the first image to be the first plurality of reference points; 
and 

step (b) includes the steps of: 

bl) identifying a second subset of features of the prese- 
lected set of features of the scene that are visible in the 
second image; and 

b2) determining depictions of the second subset of features 
in the second image to be the second plurality of refer- 
ence points. 

9. The method according to claim 7, wherein the set of 
reference points common to the first plurality of reference 
points and the second plurality of reference points includes 
exactly two, three, or four elements. 

10. The method according to claim 7, wherein: 

the first image and the second image are high altitude 
images of a planet surface; and 

the first plurality of reference points and the second plural- 
ity of reference points are ground control points (GCPs). 

11. The method according to claim 10, wherein: 

the high altitude images of the planet surface each include 
a section of coastline; and 

the GCPs are depictions of features of the section of coast- 
line. 

12 . A non-transitory computer readable medium adapted to 
instruct a general purpose computer to register a first image to 
a second image using a similarity transformation, the first 
image including a first plurality of pixels mapped to a set of 
first image coordinates and the second image including a 
second plurality of pixels mapped to a set of second image 
coordinates, the method comprising the steps of: 

a) determining the first image coordinates of a plurality of 
reference points in the first image; 

b) determining the second image coordinates of the plural- 
ity of reference points in the second image; 

c) performing a Cartesian translation of the set of second 
image coordinates such that the second image coordi- 
nates of a first reference point of the plurality of refer- 
ence points match the first image coordinates of a first 
reference point; and 

d) if the plurality of reference points includes exactly two 
reference points, performing a similarity transformation 
of the translated set of second image coordinates to scale 
the second image coordinates and rotate the second 
image coordinates about the first reference point such 
that the second image coordinates of a second reference 
point of the plurality of reference points match the first 
image coordinates of the second reference point; and 
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e) if the plurality of reference points includes more than 
two reference points and no less than two reference 
points, performing an affine transformation of the trans- 
lated set of second image coordinates to register the 
5 second image coordinates of the plurality of reference 

points to the first image coordinates of the plurality of 
reference points. 

13. The non-transitory computer readable medium accord- 
ing to claim 12, wherein step (a) of the method includes the 

to steps of: 

al) determining initial first image coordinates of the first 
reference point in the first image; 
a2) performing a Cartesian translation of the set of first 
image coordinates such that the first image coordinates 
15 of the first reference point are (0, 0); and 

a3) determining the first image coordinates of the second 
reference point in the second image using the translated 
set of first image coordinates. 

14. The non-transitory computer readable medium accord- 
20 ing to claim 12, wherein step (a) of the method includes the 

steps of: 

al) identifying a feature visible in both the first image and 
the second image to be the first reference point; 
a2) identifying another feature visible in both the first 
25 image and the second image to be the second reference 
point; and 

a3) determining the first image coordinates of the first 
reference point and the second reference point in the first 
image. 

30 15. The non-transitory computer readable medium accord- 

ing to claim 12, wherein: 

the first image is a reference image including a plurality of 
predetermined features; and 
step (a) of the method includes the steps of: 

35 al) identifying one of the plurality of predetermined fea- 
tures in the second image to be the first reference point; 
a2) identifying another of the plurality of predetermined 
features in the second image to be the second reference 
point; and 

40 a3) determining the first image coordinates of the first 

reference point and the second reference point in the first 
image. 

16. The non-transitory computer readable medium accord- 
ing to claim 12, wherein: 

45 the first image and the second image are high altitude 
images of a planet surface; and 
the first reference point is a first ground control point 
(GCP) and the second reference point is a second GCP. 

17. The non-transitory computer readable medium accord- 
50 ing to claim 16, wherein: 

the high altitude images of the planet surface each include 
a section of coastline; and 

the first GCP is a feature of the section of coastline and the 
second GCP is another feature of the section of coast- 
55 line. 



